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Abstract. The numerical studies of the interplanetary coupling between 
multiple magnetic clouds (MCs) are continued by a 2.5-dimensional ideal mag- 
netohydrodynamic (MHD) model in the heliospheric meridional plane. The 
interplanetary direct collision (DC) / oblique collision (OC) between both 
MCs results from their same/different initial propagation orientations. Here 
the OC is explored in contrast to the results of the DC (Xiong et al., 2007). 
Both the slow MCI and fast MC2 are consequently injected from the dif- 
ferent heliospheric latitudes to form a compound stream during the inter- 
planetary propagation. The MCI and MC2 undergo contrary deflections dur- 
ing the process of oblique collision. Their deflection angles of \S6i\ and \692\ 
continuously increase until both MC-driven shock fronts are merged into a 
stronger compound one. The \S9i\, \S92\, and total deflection angle A9 {A9 = 
\S9i \ + I (^^2 1) reach their corresponding maxima when the initial eruptions 
of both MCs are at an appropriate angular difference. Moreover, with the 
increase of MC2's initial speed, the OC becomes more intense, and the en- 
hancement of 69i is much more sensitive to 692- The \S9i\ is generally far less 
than the \S92\, and the unusual case of \S9i\ ~ |<56'2| only occurs for an ex- 
tremely violent OC. But because of the elasticity of the MC body to buffer 
the collision, this deflection would gradually approach an asymptotic degree. 
As a result, the opposite deflection between the two MCs, together with the 
inherent magnetic elasticity of each MC, could efficiently relieve the exter- 
nal compression for the OC in the interplanetary space. Such deflection ef- 
fect for the OC case is essentially absent for the DC case. Therefore, besides 
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the magnetic elasticity, magnetic helicity, and reciprocal compression, the 
deflection due to the OC should be considered for the evolution and ensu- 
ing geoeffectiveness of interplanetary interaction among successive coronal 
mass ejections (CMEs). 
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1. Introduction 

One of the greatest concerns within the current space science community has been 
increasingly focused on the Sun-Earth system, which is intimately linked by the solar 
wind. The solar wind originates from the chromospheric network [Xia et ai, 2003; Xia, 
2003; Xia et ai, 2004], according to the measurements of ultraviolet emission and Doppler 
shifting speed in the inner corona, carries non-Wentzel-Kramers-Brillouin Alfven Waves 
in the differential flow of multiple ion species [Li and Li, 2007a, 2008], is very likely driven 
by an ion cyclotron resonance mechanism via the Kolmogorov turbulent cascade [Li et ai, 
2004], and transports the angular momentum from the Sun [Li and Li, 2006; Li et ai, 
2007b; Li and Li, 2009]. The ubiquitous interplanetary solar wind highly fluctuates, owing 
to outward-emanating disturbances from solar activities. Therefore, the Sun serves as the 
driver for the cause-and-effect transporting chain of space weather. 

The interplanetary space, which Dryer [1994] calls a "transmission channel" between 
the Sun and the Earth , is a nonlinear system consisting of various discontinuous fronts, 
diffusion processes, and couplings between different spatial and temporal scales. A mag- 
netic flux rope levitating in the corona may suddenly lose its equilibrium and consequently 
escape into the interplanetary space [Chen et ai, 2006, 2007]. The interplanetary mani- 
festation of such a magnetic rope is identified as a magnetic cloud (MC) with enhanced 
magnetic field magnitude, smooth rotation of the magnetic field vector, and low proton 
temperature [Burlaga et ah, 1981]. The passage of an MC across the Earth triggers a 
significant geomagnetic storm because of large southward magnetic flux in the MC body 
[Tsurutani et ai, 1988; Gosling et ai, 1991]. Hence MCs are an important subset of 
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interplanetary CMEs (ICMEs), whose fraction is ~ 100%, though with low statistics, at 
solar minimum and ~ 15% at solar maximum [Richardson and Cane, 2004]. Especially 
at solar maximum, when the daily occurrence rate of CMEs is about 4.3 on average based 
on the SOHO/Lasco CME catalogue (http://cdaw.gsfc.nasa.gov/CME_list), CMEs very 
likely interact with each other on their journey toward the Earth. Two distinct obser- 
vation events of interaction between an early slow CMEl and a late fast CME2 within 
30 solar radii were presented by Gopalswamy [2002]: (1) two fast CMEs on 4 Novem- 
ber 1997, which were initially 100° apart in relation to their source regions according 
to Yohkoh/SXT observation, led to the plowing of the CME2-driven shock through the 
CMEl in the field-of-view (FOV) of Lasco C2/C3; (2) two fast CMEs from the same 
source region on 20 January 2001, initially two hours apart, were later indistinguishable 
in the FOV of Lasco C3, and are therefore thought to have cannibalized each other. As 
the coupling of multiple CMEs from the same/different heliographic location of source 
region is defined as the direct collision (DC) / oblique collision (OC) by Xiong et al. 
[2006b], these two events of 4 November 1997 and 20 January 2001 are the cases of OC 
and DC, respectively. The radio signatures of coronal mass ejection cannibalism typi- 
cally precede the intersection of the leading-edge trajectories and behave as an intense 
continuum-like radio emission enhancement, usually following a type II radio burst on 
basis of Wind / WAV E S ohsQi\SitioYi [Gopalswamy et al., 2001, 2002]. Gopalswamy et al. 
[2002] argue that the nonthermal electrons responsible for this new type of ratio emission 
are accelerated due to magnetic reconnection between two CMEs and/or the formation of 
a new shock at the time of collision between two CMEs. Meanwhile, some interplanetary 
complicated structures were also reported in the near-Earth space, such as the complex 
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ejecta [Burlaga et ai, 2002], compound stream [Burlaga et ai, 1987; Wang et ai, 2003a; 
Dasso et ai, 2009], shock-penetrated MCs [Lepping et ai, 1997; Wang et ai, 2003b; 
Berdichevsky et ai, 2005], and non-pressure-balanced "MC boundary layers" associated 
with magnetic reconnection [Wei et ai, 2003, 2006]. According to the in situ observations 
of spacecraft at 1 AU, the evolutionary signatures of ICMEs' interaction include heating 
of the plasma, acceleration/deceleration of the leading/trailing ejecta, compressed field 
and plasma in the leading ejecta, possible disappearance of shocks, and strengthening 
of the shock driven by the accelerated ejecta [Farrugia and Berdichevsky, 2004]. Since 
magnetic diffusion in interplanetary space is much less than that in the solar corona, the 
cannibalism of CMEs that interact in the Lasco FOV [Gopalswamy et ai, 2001, 2002] 
should not occur in the interplanetary space [Xiong et ai, 2007]. Moreover, formed by 
multiple CMEs/ICMEs colliding, the compound stream at 1 AU could be in a different 
evolutionary stage. The position of the overtaking shock at 1 AU can be (1) still in 
the MC, such as an 18 October 1995 event [Lepping et ai, 1997] and a 5-7 November 
2001 event [Wang et ai, 2003b], or (2) ahead of the MC after ultimately penetrating 
it [Berdichevsky et ai, 2005]. The compressed magnetic field downstream of the shock 
front is northward for the 18 October 1995 event [Lepping et ai, 1997] and southward 
for the 5-7 November 2001 event [Wang et ai, 2003b]. Therefore, the latter event of 5-7 
November 2001 resulted in a great magnetic storm of Dst ~ —300 nT. An important 
interplanetary origin for the great geomagnetic storms have already been identified by the 
observations [Wang et ai, 2003a; Farrugia et ai, 2006; Dasso et ai, 2009] and simula- 
tions [Xiong et ai, 2006a, b, 2007; Xiong, 2007] as multi-ICME structures, accompanying 
intense compression of southward magnetic flux during the interaction process. When 
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the compound structure reaches the Earth through the interplanetary space, its physi- 
cal parameters are jointly decided by three factors: (1) individual CMEs themselves, (2) 
inhomogeneous interplanetary medium, (3) irreversible interacting process among these 
CMEs/ICMEs [Xiong, 2007]. Due to the intractability of analytical reduction, compound 
structures resulting from the interaction of multiple CMEs/ICMEs have been extensively 
studied in numerical simulations: e.g., complex ejecta [Xiong et al, 2005], interaction of 
a shock wave with an MC [Vandas et al, 1997; Xiong et ai, 2006a, b], and coupling of 
multiple MCs [Schmidt and Cargill, 2004; Lugaz et a/., 2005; Xiong et a/., 2007]. Partic- 
ularly, Xiong et al. [2005, 2006a, b, 2007] and Xiong [2007] conducted a systematic and 
delicate numerical MHD simulation of interplanetary compound structures in terms of 
their formation, propagation, evolution, and ensuing geoeffectiveness. These simulation 
works do well provide theoretical interpretations for physical phenomena of compound 
structures observed by the SOHO, Wind, and ^4 Ci? spacecraft. 

The radial lift-off of a CME at its onset phase from a solar source region sometimes 
deviates from the radial ray during its outward movement. The non-straight trajectory 
substantiates that deflections do happen during CME/ICME propagation. The deflection 
effect plays a notable role in space weather predicting, since the first step of predic- 
tion is whether or not a solar eruption will ultimately affect the geospace environment 
[Williamson et al, 2001]. The near-Sun trajectory of a CME can be directly imaged by 
remote sensing of a white light coronagraph onboard such spacecraft as Skylab, SOHO, 
and STEREO. MacQueen et al. [1986] found that 29 CME events observed during the 
Skylab epoch of solar minimum from 1973 to 1974 underwent an average 2.2° equator- 
ward deflection, and ascribed that the deflection to the nonradial forces arising from the 
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background coronal magnetic and flow patterns. Cremades and Bothmer [2004] identified 
the CME events from SOHO/Lasco FOV and their corresponding source regions from the 
SOHO/EIT and SOHO/MDI from January 1996 to December 2002, and found that the 
position angle (PA) of Lasco-imaged CMEs deviates statistically about 18.6° southward 
toward the lower latitude at solar minimum. They also ascribed such equatorward de- 
flection of CMEs from solar activity belts to the surrounding fast solar wind from polar 
coronal holes with a stronger total plasma and magnetic field pressure. Gopalswamy et al. 
[2001] reported that on 10 June 2000, a slow CME of 290 km/s was overtaken by a fast 
CME of 660 km/s from a different solar source region; the core of the slow CME was 
leftward deviated by 13° in terms of the PA in the Lasco/C3 FOV. Zhang et al. [2004] 
also reported a nonradial motion of a gradually accelerated CME on 19 October 1997 from 
the SOHO observation. This peculiar CME was initiated above the east limb at northern 
latitude 14°N in the EIT FOV, tilted towards the equator as it rose in the Lasco/Cl 
FOV, and was very symmetric with respect to the equator later in the Lasco C2/C3 
FOV. Furthermore, besides the occurrence within the Lasco/C3 FOV, the CME/ICME 
deflection does exist beyond the near-Sun space. On the basis of statistical analyses of 
interplanetary scintillation observations, Wei [1988] and Wei and Dryer [1991] found that 
the solar-flare-generated shock deflects eastward in the heliospheric equator and equator- 
ward in the heliospheric meridian during its interplanetary propagation. This deflection 
evidence of interplanetary shock aphelion results from joint effects of the (1) spiral in- 
terplanetary magnetic field (IMF), (2) westward movement of the heliographic location 
of a solar flare during the impulsive phase, and (3) heterogenous medium consisting of 
the fast solar wind from open corona magnetic field and the slow solar wind astride the 
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heliospheric current sheet [Hu, 1998; Hu and Jia, 2001]. In addition, the solar source 
distribution of Earth-encountered halo CMEs is east- west asymmetry [Wang et ai, 2002; 
Zhang et ai, 2003]. Some eastern limb CMEs hit the Earth [Zhang et ai, 2003], and 
conversely some disk CMEs missed the Earth [Schwenn et ai, 2005]. According to an 
ICME's kinematic model [Wang et ai, 2004], ICMEs could be deflected as much as sev- 
eral tens of degrees during its propagation by the background solar wind and spiral IMF; 
a fast CME will be blocked by the background solar wind ahead and deflected to the east; 
a slow CME will be pushed by the following background solar wind and deviated to the 
west. The existence of ICME deflection is obviously implied from the evidence of indirect 
observations about the correlations between the near-Sun CME and near-Earth ICME. 
However, direct observations covering the entire interplanetary space have only been avail- 
able since the launching of SMEI and STEREO in the twenty-flrst century. Most of the 
current spaceborne observations are still heavily concentrated to the thirty solar radii by 
remote sensing, and the geospace by in situ detecting. As interplanetary observation data 
is relatively small, numerical simulations are necessary and signiflcant for understanding 
the whole of interplanetary dynamics, including the deflection effect. Xiong et al. [2006b] 
proposed that (1) the OC between a preceding MC and a following shock results in the 
simultaneous opposite deflections of the MC body and shock aphelion; (2) an appropriate 
angular difference between the initial eruption of an MC and an overtaking shock leads 
to the maximum deflection of the MC body; (3) the larger the shock intensity is, the 
greater the deflection angle. As a straightforward analogy to the MC-shock OC [Xiong 
et al., 2006b], the interplanetary deflection can be also expected for the MC-MC OC. As 
a result of collision of one MC with either a shock or another MC, the deflections can be 
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ascribed to the interaction between different interplanetary disturbances. In contrast to 
our models, the previous deffection models [e.g., Hu, 1998; Hu and Jia, 2001; Wang et ai, 
2004] are caused by the interaction between the ambient solar wind and interplanetary 
disturbance. 

The conjecture about the interplanetary deffection from the MC-MC OC is investigated 
in this paper. In addition, a simplified circumstance of MC-MC DC, excluding the deflec- 
tion effect, has already been studied by Xiong et al. [2007]. The following conclusions are 
revealed from the MC-MC DC [Xiong et al, 2007]: (1) when the accumulated magnetic 
elasticity can balance the external colliding, the compressibility of double MCs reaches its 
maximum; (2) this cutoff limit of compressibility mainly decides the maximally available 
geoeffectiveness of double MCs, because geoeffectiveness enhancement of MCs' interacting 
is ascribed to compression; (3) the magnetic elasticity, magnetic helicity of each MC, and 
compression between each are the key physical factors for the formation, propagation, 
evolution, and resulting geoeffectiveness of interplanetary double MCs. Here the study of 
MC-MC OC is a more reasonable extension of that of MC-MC DC [Xiong et al, 2007]. 
Thus two issues are naturally raised: (1) What is the difference between the MC-MC DC 
and MC-MC OC in terms of the interplanetary dynamics and ensuing geoeffectiveness? 
(2) Does such a deflection effect caused by the MC-MC OC play a significant or negli- 
gible role during interaction process? The answers to these questions are explored by a 
2. 5- dimensional (2.5-D) numerical model in ideal MHD process. 

The present work targets the OC between two MCs as our logical continuation in a series 
of studies for the interplanetary compound structures [Xiong et al., 2005, 2006a, b, 2007; 
Xiong, 2007]. We give the numerical MHD model in section 2, describe the dynamics 
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and geoeffectiveness of two typical cases of double MCs in section 3, analyze the roles of 
eruption interval in section 4, angular difference in section 5, and collision intensity in 
section 6 for two MCs' interacting, and summarize the paper in section 7. 

2. Numerical MHD Model 

The dynamics and geoeffectiveness of interplanetary compound structures have al- 
ready been numerically investigated by our effective numerical model [Xiong et ai, 
2006a, b, 2007; Xiong, 2007]. This model quantitatively relates the output of solar dis- 
turbances at 25 Rg to the interplanetary parameters and geomagnetic storm at 1 AU, 
thus establishing a cause-and-effect transporting chain for a solar-terrestrial physical pro- 
cess. The concrete implementation of this numerical model consists of two steps: (1) the 
numerical MHD simulation of interplanetary disturbance propagation, and (2) using the 
Burton empirical formula for the solar wind - magnetosphere - ionosphere coupling to 
evaluate the geomagnetic storm index Dst [Burton et ai, 1975]. The detailed description 
of the numerical model, including the numerical algorithm, computational grid layout, 
ambient solar wind, is given in Xiong et al. [2006a]. 

An incidental MC, radially launched from the solar surface, is characterized by several 
parameters: the emergence speed Vmc, latitude 9mc, and time tmc, et al. The following 
MC2's emergence latitude 9mc2 is included for parametric study in contrast to the DC along 
the heliospheric equator [Xiong et al, 2007]. Both the MCs are consequently injected into 
the simulation domain through a particular modification of the inner boundary condition 
at 25 Rs [Vandas et ai, 1995; Xiong et ai, 2006a]. The DC and OC in the interplanetary 
medium correspond to 9m.c2 = 0° and 9mc2 7^ 0°, since the preceding MCI emerges from 
the equator, 9mci = 0°. Moreover, the MC2-driven shock in all of our simulation cases 
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is faster than the local magnetosonic speed at any time in order to prevent weak shock 
dissipation in the MC body of low plasma (3. 

3. MC1-MC2 Interaction 

All thirty-four cases of double MCs' interacting are assembled into three groups in Table 
2, with ten cases of an individual MC in two groups from Table 1 for comparison. In Table 
2, the helicity of one MC i/mci = 1 is opposite to that of the other MC Hmc2 = because 
we are only interested in the maximally available geoeffectiveness among all combinations 
of each MC helicity [Xiong et ai, 2007]. Groups of an individual-preceding MC (IPM), 
an individual-following MC (IFM), an eruption- interval dependence (EID), an angular- 
difference dependence (ADD), and a collision-intensity dependence (CID) are studied, 
with case E2 shared by Groups EID, ADD, and CID. The slow MCI of Vmci = 400 
km/s, Hmci = 1, ^mci = 0°, and tmci = hour is chased and pounded by a fast MC2 of 
various parameters. Here the parametric studies of double MCs cover a wide spectrum of 
tmc2 = 10.2 ~ 44.1 hours in group EID, 6^02 = 0° ~ 50° in group ADD, and Vmc2 = 450 ~ 
1200 km/s in group CID. Moreover, by adjusting Dt {Dt = t^ci — tmc2, tmci = hour), 
the initiation delay between the two MC emergences in group EID, an interplanetary 
compound stream consisting of double MCs may reach a different evolutionary stage 
when it arrives at 1 AU. The tmc2 is prescribed to be 12.2 hours in groups ADD and 
CID for the full development of double MCs' interacting within 1 AU. In the following we 
address case Ei of 30.1 hours and case E2 of 12.2 hours in group EID, which are typical 
examples of double MCs in the early and late evolutionary stages. 

3.1. Case Ei 
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Figures 1-4 shows the consequent behavior of MC1-MC2 interaction of Case Ei with 
the eruption speed Vmci = 400 km/s, Vmc2 = 600 km/s, and the initiation delay tmc2 = 
30.1 hours. The magnetic field lines, of which two are enclosed with white solid lines 
marking the boundaries of MCI and MC2, are superimposed on each image of Figures 
1-3. Two radial profiles, one through the equator (noted by Lat. = 0°), the other through 
4.5° southward (white dashed lines in Figures 1-3, noted by Lat. = 4.5°S), are plotted 
in Figure 4. The magnitude B of the magnetic field in the radial profile of Figures la-c 
is presented by subtracting its initial value B\t=o of ambient equilibrium. The coupling 
of two MCs could be considered a comprehensive interaction between two systems, each 
comprised of an MC body and its driven shock. The MC2-driven shock and MC2 body 
are successively involved in the interaction with the MCI body. The MC2-driven shock 
catches up with the MCI body tail at 48 hours, as seen in Figures 2d and 4d. Across 
the shock front, impending collision is influenced by the abrupt jump of radial speed Vr 
from 430 to 650 km/s. From then on, both MCs are coupled with each other to form an 
interplanetary compound stream of double MCs [Wang et al, 2003a; Dasso et al, 2009]. 
At 57 hours, the marching MC2-driven shock front behaves as a steep speed jump at MCl's 
rear part (Figures 2e and 4e), just downstream from which the magnetic magnitude B 
(Figures lb and 4b) and fast magnetosonic mode speed Cf (Figures 3h and 4h) are locally 
enhanced. Due to the large initial delay tmc2 = 30.1 hours, only the rear half of the MCI 
body is swept and compressed by the MC2-driven shock within 1 AU (Figures 2f and 4f). 

The in situ observation along Lat. = 4.5°S by a hypothetical spacecraft at the La- 
grangian point (LI) is shown in Figure 5. The boundary and core of each MC are identi- 
fied as dashed and solid lines, respectively. The rear half of the MCI body is significantly 
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gripped by the penetration of the MC2-driven shock at the MCI core and the push of the 
MC2 body upon the MCI tail. The duration of MCl's rear half (9 hours) is much less 
than that of MCl's anterior half (16 hours). The dawn-dusk electric field VB^ is swiftly 
intensified from at 73 hours to —13 mV/m at 76 hours. Because the orientation of 
the magnetic field within the double flux-rope structure is north-south-south-north, the 
superposition of three individual southward Bg regions from the MCI, IMF, and MC2 
behaves as a long-lived geoeffective solar wind flow from 73 to 93 hours (Figure 5d), and 
results in a one-dip curve of Dst with its minimum —234 nT at 87 hours (Figure 5e). 

3.2. Case E2 

In case E2, a much earlier emergence time of the MC2 (tmc2 = 12.2 hours) guarantees 
the full interaction between the two MCs before their arrival at 1 AU. Only the evolution 
of Vr is given in Figures 6 and 7 to visualize the structure of double MCs. The initial 
emergence latitudes of MCI {Omci\t=o) and MC2 {6mc2\t=o) are two important parameters 
of solar eruption output. The nonzero difference D9\t=o {D6\t=o = 9mc2\t=o — ^»nci|t=o) 
decides a consequent OC in the interplanetary space. The collision between two MCs 
can be understood by comparison to a billiards game. For one moving rigid ball colliding 
with another still ball along the radial direction, the response is straightforward in a 
vacuum: in the DC case, two balls will strictly move along the same radial direction; in 
the OC case, two balls will oppositely deflect along an angular direction, accompanying 
their continuous radial movement. The patterns of ball movement contribute to a further 
quantitative understanding of complex collision between two MCs in the interplanetary 
medium. The magnetic field lines, frozen in a low P plasma, could be considered as an 
elastic skeleton embedded in the MC body. The innate magnetic elasticity can efficiently 
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buffer the compression as a result of the colhding of a following MC2 against a leading MCI 
[Xiong et ai, 2007]. When every MC becomes increasingly stiff, the compression reaches 
its asymptotic degree. The compressibility effect should be included for the quantitative 
investigation of deflection effect as a result of OC. Such a task should use numerical 
MHD simulation, as we demonstrate in this paper. The direction of main compression 
within the double MCs is parallel/oblique to the radial direction for the DC/OC. For the 
DC already discussed by Xiong et al. [2007], the compression strictly persists along the 
heliospheric equator and the compressed magnetic flux within the MC body almost points 
to the south. Therefore, the DC case is very efficient to enhance the geoeffectiveness. The 
fast MC2 body continuously strikes the slow MCI tail, until the MC2 speed is lower than 
the MCI speed after momentum transfer. Such MC2 body pushing prevents magnetic 
field lines in the MCI, previously compressed by MC2-driven shock, from being restored, 
when the MC2-driven shock completely passes through the MCI body. For the OC, the 
compression occurs along one side of each MC. For an example, the MC2 body is faced 
with the MCI body from its left side and the ambient solar wind from its right side. 
Due to the MCl's blocking, the MC2 suffers the compression from its left side. Such 
an angular pressure imbalance leads to the MC2's right deflection. Simultaneously, the 
MCI deflects leftward for the same reason. The opposite deflections, separating double 
MCs, greatly relieve the intensity of the OC. Therefore, the angular freedom for each 
MC is an extra factor in efficiently buffering the compression. This angular deflection, 
absent for the DC case [Xiong et al., 2007], is explored for the OC case here. The IMF 
lines within the latitude difference {D6\t=Q = 10° in this case) of two MC eruptions are 
first draped and then compressed between the MCI tail and MC2 head. At 22 hours. 
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the left flank of the MC2-driven shock enters the MCI core (Figures 6a and 7a). Due 
to very low /9 in the MCI medium, the left flank of the shock front in the MCI body 
propagates much faster than the right one in the ambient solar wind. The MCI body is 
compressed by the MC2-driven shock along its normal. The advance of the MC2-driven 
shock accompanies a drastic jump of local speed in the MCI medium. The MC2 body 
obliquely chases the MCI body and then grazes the MCl's right boundary. During this 
process, the momentum is gradually transferred from the following MC2 to the preceding 
MCI. The location of the most violent interaction within the double MCs, characterized 
by the greatest compression of local magnetic flux, gradually shifts from the MCl's rear 
half (Figures 6a and 7a) to the MCl's anterior half (Figures 6b and 7b), and is finally 
within the MCl-driven sheath (Figures 6c and 7c). The compound stream of double MCs 
reaches a relatively stable state at 57 hours (Figures 6c and 7c) when the MC2-driven 
shock ultimately merges with the MCl-driven shock into a stronger compound one. 

The time sequence of synthetic measurement at LI for case E2 is shown in Figure 8. 
The speed Vr monotonically decreases from the MCl's head to the MC2's tail (Figure 
8c). The magnetic elasticity of southward magnetic flux takes a recovering effect against 
the previous compression, as the MC2-driven shock continuously moves forward in the 
MCI body. As the compression of the MCl's rear half is largely relieved, the duration of 
geoeffective solar wind flow is prolonged from 20 hours in case Ei to 31 hours in case E2. 
Owing to the OC, the MCI deflects northward, and the MC2 deviates southward. Largely 
reduced is the total southward magnetic flux passing through Lat. = 4.5°S. The opposite 
deflection of the two MCs together with the above-mentioned mitigated compression cause 
a significant increase of Dst from —234 to —121 nT. 
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The evolution of various physical parameters for each MC in case E2 is shown in Figure 
9. The MCI is accelerated and the MC2 is decelerated, as seen in Figure 9a. The 
MCI begins to deflect northward at 28 hours, 16 hours later than the MC2's southward 
deflection (Figure 9b). The deflections of both MCI and MC2 gradually approach an 
asymptotic values SOmci = —3° and 66mc2 = 6.3°, respectively. Last, after being pushed 
aside, the MCI and MC2 propagate along the latitudes of 9mci = 3°N and 9mc2 = 16.3°S, 
respectively. Obviously, the MC2 undergoes a larger deflection than the MCI. Due to 
the deflection, the distance d between the two MC cores is highly increased in contrast to 
an uncoupling case (Figure 9c). For an individual MC, the higher the MC speed is, the 
greater the compression between the MC body and its front ambient solar wind, and the 
smaller the MC's cross-section area Amc- For an OC case of double MCs, Amc depends on 
one more compression factor, interaction between two MCs. With the increased speed, 
the MCI suffers larger compression from its front ambient solar wind. The MC2's trailing 
pounding the MCI compresses the MCI body. Therefore, the MCI area ^mci is smaller 
than its corresponding isolated case (Figure 9d), which is consistent with the DC case 
[Xiong et ai, 2007]. However, the MC2 area Amc2 for the OC (Figure 9e) is quite contrary 
to that for the DC [cf. Figure 5d in Xiong et ai, 2007]. The inconsistency is ascribed 
to the deflection in the OC. For the DC, the compression of MC2 chiefly exists between 
the MCI tail and the MC2 head. The persistent blocking of the MCI body causes the 
MC2 area shrinkage. For the OC, the MC2 senses the compression from two aspects: 
(1) the front solar wind, and (2) the sideward MCI body. The MC2 slowdown tends to 
enhance Amc2, which indicates the mitigated compression between the MC2 body and 
the solar wind. The blocking of the MCI body at MC2's left tends to reduce Amc2- As 
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the MC2 deflects sideward, the former factor of the front solar wind contributes more to 
Ajnc2i and the latter factor of the sideward MCI body contributes less. The integration 
of both competing factors determines the increase or decrease of Amc2 in contrast to its 
corresponding individual MC case. As for the current case E2, Amc2 is increased. 

3.3. Latitudinal Distribution of Geoeffectiveness 

It has been substantiated from both observation data analyses [e.g., Burlaga et ai, 
1987; Wang et ai, 2003a; Farrugia et ai, 2006; Dasso et ai, 2009] and numerical simu- 
lations [e.g., Xiong et ai, 2006a, b, 2007; Xiong, 2007] that multiple ICME interactions 
can significantly enhance the geoeffectiveness at 1 AU. The near-equator latitudinal dis- 
tribution of Dst index is plotted in Figure 10. Since an individual MC would propagate 
radially through interplanetary space, the MC core passage corresponds to the strongest 
geomagnetic storm in a one-dip latitudinal distribution of geoeffectiveness. The strongest 
geoeffectiveness is —103 nT at 0° for an isolated MCI event (case Pi) and —140 nT at 
10°S for an isolated MC2 event (case F3). The coupling of two MCs obviously aggravates 
the geoeffectiveness. For case Ei, the geoeffectiveness of the two MCs is overlapped, so 
that the 9 — Dst curve looks like a single dip with its minimum —262 nT at 9°S. For case 
E2, the initial delay between the two MCs is short {tmc2 = 12.2 hours), the double MCs ex- 
perience sufficient evolution, the accumulated deflection angle becomes very pronounced, 
the latitudinal distance between the two MCs becomes large, and the geoeffectiveness of 
the two MCs is thus separated; hence the 6 — Dst curve behaves like two local dips with 
their local minima of —200 nT at 15°S and —145 nT at 1.5°N. As the compound stream 
at 1 AU formed by the two MCs' coupling evolves from case Ei to E2, its geoeffectiveness 
is significantly diffused along the latitude with the intensity largely reduced. 
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The interplanetary dynamics and resulting geoeffectiveness of the double MCs is a com- 
plex system involving multiple independent variables. The parametric studies of eruption- 
interval dependence (EID), angular-difference dependence (ADD), and collision-intensity 
dependence (CID) are further explored below to continue our preliminary efforts for the 
cases El, E2 of double MCs, and the cases Pi, F3 of a single MC. In sections 4, 5, and 
6, the geoeffectiveness of double MCs is described by a scalar of minimum Dst along its 
latitudinal distribution. 

4. Eruption-Interval Dependence 

The results of eruption- interval dependence is elucidated in Figure 11. The transfer 
of momentum from the fast MC2 to the slow MCI leads to shortening of the Sun-Earth 
transient time TTmci and lengthening of TT„^c2 (Figure 11a). As the initial eruption delay 
{tmc2 — tmci) IS shortened, the deflection of each MC exhibits an asymptotic behavior 
(Figure lib). When tmc2 is reduced from 33.1 to 22.1 hours, and then to 10.2 hours, the 
MCI deflection angle |5^mci| is increased from 0.2° to 0.7°, and then to 3.3°; the MC2 
deflection angle |5^mc2| is enhanced from 1.2° to 2.7°, and then to 6.9°; the total deflection 
angle AO {AO = \60mci \ + \S6mc2\) is changed from 1.4° to 3.4°, and then to 10.2°. These 
deflection ratios of |5^mci| : |5^mc2| at tmc2 = 33.1, 22.1, and 10.2 hours correspond to 
0.17, 0.26, and 0.48, respectively. Obviously, the MC2 occupies a much bigger share of 
the total deflection angle A6. This latitudinal deflection is manifested in the distance 
d between the cores of MCI and MC2 at 1 AU (Figure lid). At the beginning, the 
magnetic field lines in the MCI rear half are too vulnerable to resist the MC2's pounding, 
so the temporarily enhanced compression leads to d decrease. As interpreted in section 
3.2, when the most intensely interacting region in the double MCs is shifted from the 
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MCl's rear half, the magnetic elasticity, being passively quenched at the earlier time by 
the MC2-driven shock, begins to actively bounce to push the two MCs apart. Then the 
d is steadily increased. Hence the minimum d of 63 Rs exists in an intermediate platform 
of tmc2 = 22 ~ 28 hours. Between this zone of tmc2 = 22 ~ 28 hours, the double MCs 
suffer the strongest compression, the MCl's cross section area Amci is compressed to its 
minimum 2.4 x 10'^ i?^ at tmc2 = 22 hours (Figure lie), and Dst reaches a minimum of 
—270 nT at tmc2 = 28 hours (Figure llf). When tmc2 < 22 hours, the magnetic elasticity 
restoration and angular deflection lead to the increasingly weak compression between 
two clouds and a consequent increase in Amc\- In addition, according to the reasons 
given in section 3.2 for the Amc2 variance, the predominance of the MC2's momentum 
loss over the MCl's blocking is responsible for the monotonic increase of Amc2 between 
Dt = —32 ~ — 10 hours. These behaviors of the OC case are essentially different from 
those of the DC case [Xiong et ai, 2007]. For the DC case without the deflection effect, 
the persistent following of the MC2 body at the MCI tail can inhibit the MCI body from 
re-expanding, so the Dst can be roughly maintained at a constant, given that the initial 
delay tmc2 is smaller than a certain threshold [cf. Figure 9 in Xiong et ai, 2007]. The 
perpetual balance between the external compression and innate elasticity for the DC case 
is out of equilibrium for the OC case under a new circumstance of angular deflection. 
Moreover, the external compression could be largely offset by the deflection. Therefore, 
the OC case is generally weaker for geoeffectiveness than its corresponding DC case, and 
the strongest geoeffectiveness for the OC case can only be achieved at a certain initial 
delay between two MC eruptions. 
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5. Angular-Difference Dependence 

The angular-difference dependence is shown in Figure 12. The initial latitudinal dif- 
ference of D6\t=Q = Omc2\t=o — (^mci\t=o for solar output decides the oblique degree of 
interplanetary collision between two MCs. Corresponding to the nonexistence of deflec- 
tion, D9\t=Q = 0° is due to the symmetrical condition, which was thoroughly addressed 
by Xiong et al. [2007]. When D9\t=Q is too large, the OC effect will be significantly 
mitigated, and the consequent deflection will be obviously weak. An appropriate D6\t=Q 
corresponds to the maximum deflection of OC cases of double MCs, very similar to the 
known conclusion for the OC case of "a shock overtaking an MC" [Xiong et al., 2006b]. At 
D9\t=Q = 15°, the total deflection angle A6 reaches its maximum 12.2° with 66mci = —3.8° 
and 69mc2 = 8.4°. The |(56'mc2| is generally larger than the |(56'mci|, but it does not match 
the case of D9\t=o > 40°. When D9\t=o > 40°, the two MCs are so widely separated that 
the interaction is virtually ascribed to the coupling of the MCI body and the MC2-driven 
shock. Such indirect interaction between the two MC bodies to transfer momentum clar- 
ifles l^^mcil > |f^6'mc2| for D9\t=o = 40° ~ 50° (Figure 12b) and the Amc2 decrease for 
D9\t=o = 20° ~ 50° (Figure 12e). With respect to the Amc2 variance for D9\t=o < 20°, 
two competing factors of MC2 momentum loss and MCI body blocking take effect, as pre- 
viously interpreted in section 3.2. The dominance of the MC2 momentum loss accounts 
for the increase between D9\t=o = 10° ~ 20°; that of the MCI body blocking elucidates 
the decrease between D9\t=o = 0° ~ 10°. The closer the two MCs are in the near-Sun 
position, the smaller the distance d at 1 AU (Figure 12d). As -06*14=0 decreases, the Dst, 
being steadily reduced with a steeper slope, changes from —140 nT at D9\t=o = 50° to 
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—230 iiT at D9\t=o = 0°. The less the deflection effect is, the more compact the multiple 
interplanetary geoeffective triggers and the more violent the ensuing geomagnetic storm. 

6. Collision-Intensity Dependence 

Figure 13 displays the collision-intensity dependence. The variance of Vmc2 corresponds 
to a different individual MC2 event. As Vmc2 increases, both TTmci and TTmc2 decrease. 
However, the decreased TTmc2 in the case of double MCs is still larger than its corre- 
sponding individual MC case (Figure 13a). The influence of the OC intensity can be 
described by the Vmc2 in some senses. As an asymptotic response to the Vmc2 increase 
from 450 to 1200 km/s, the geoeffective Dst decreases from —145 to —255 nT (Figure 
13f), and the total deflection angle increases from 7.5° to 12.8° (Figure 13c). The 
contribution of almost stems from the MCI deflection Smcii since the MC2 deflection 
angle 5mc2 is nearly constant at 6° (Figure 13b). The deflection ratio between the MCI 
and MC2 |5^mci| : |^^mc2| is 0.3 at Vmc2 = 450 km/s, 0.7 at Vmc2 = 800 km/s, and 1 at 
Vmc2 = 1200 km/s. Therefore, the cause of intensity aggravation of two MCs' colliding is 
mainly manifested in the response of the preceding MCI body. In addition, for the case of 
D9\ t=o = 10° in group ADD, the MC2's cross-section area Amc2 of a coupled case is larger 
than that of an isolate case (Figure 12e). Since D9\t=o equals 10° in group CID, the Amc2 
behavior is similar for the reason explained in section 5. Furthermore, the more intense 
the OC is between two MCs, the more violent the compression in the double MCs, and 
the stronger the accumulated innate magnetic elasticity against the external compression. 
So the deflection angle 69 (Figure 13b), the distance between the two MC cores d (Figure 
13d), and the geoeffective Dst (Figure 13f) all exhibit an asymptotic behavior. 
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7. Conclusions and Summary 

The dynamics and geoeffectiveness of interplanetary compound structures such as the 
complex ejecta [Xiong et ai, 2005], MC-shock [Xiong et ai, 2006a, b], and MC-MC [Xiong 
et ai, 2007] have been comprehensively investigated during the recent years with our 2.5- 
D numerical model within an ideal MHD framework. As a logically direct continuation to 
the DC mode between a preceding MCI and a following MC2 [Xiong et ai, 2007], the OC 
mode is further explored here to highlight a deflection effect from the parametric stud- 
ies of eruption-interval dependence, angular-difference dependence, and collision-intensity 
dependence. The deflection angle for an MC1-MC2 OC in this paper is obviously greater 
than that for an MC-shock OC addressed by Xiong et al. [2006b], as the MC1-MC2 cou- 
pling involves a comprehensive interaction among the MCl-driven shock, the MCI body, 
the MC2-driven shock, and the MC2 body 

An interplanetary compound stream is formed as a result of interaction between two 
MCs in the interplanetary space. The direction of main compression within the double 
MCs is parallel/oblique to the radial direction for the DC/OC. The OC leads to first 
compress each MC on one side, then push the MC to the other side as a result of angular 
pressure imbalance. Such a deflection effect for the OC case is essentially absent for the 
DC case. The deflection angles of MCI (|5^i|) and MC2 (|5^2|) asymptotically approach 
their corresponding limits, when the two MC-driven shocks are merged into a stronger 
compound shock. During this process, the geoeffectiveness of double MCs is significantly 
diffused along the latitudinal distribution, with the intensity largely reduced. An appro- 
priate angular difference between the initial eruptions of two MCs leads to the maximum 
deflection of \S6i\ and \662\. A continuous increase of OC intensity can synchronously 
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enhance \S6i\ and \S92\, although its effect becomes less and less obvious. The response 
of \S6i\ is far more sensitive than that of \S92\. The \S9i\ is generally far less than the 
\S92\, and the unusual case of \S6i\ ~ \S92\ only occurs for the extremely intense OC. The 
opposite deflection between two MCs, together with the inherent magnetic elasticity of 
each MC, could efficiently buffer the external compression for the interplanetary OC. 

The axial variance of an MC is ignored in our model for simplification, so that the 
geometry of an MC is reduced to be 2.5-D. In reality, both feet of an interplanetary MC 
is still connected to the solar surface, as substantiated from the evidence of bi-directional 
electron fluxes along an MCs axis [Larson et at, 1997]. However, for the local analyses of 
a cross section of an MC, a locally cylindric flux-rope has widely been used to approximate 
the globally curved one, such as the data inversion from the near-Earth in-situ observations 
[Burlaga et al., 1981], the kinematic model of an MC propagation [Owens et ai, 2006], 
and the numerical simulation of magnetic-flux-rope dynamics [Schmidt and Cargill, 2004; 
Xiong et ai, 2006a]. Hence, our 2.5-D model can well reflect some dynamic characteristics 
of 3-D MCs to some extent. 

An assimilatively integrated study of observation data analyses and numerical simula- 
tions is crucial and effective for an in-depth and overall understanding of the Sun-Earth 
system. As a CME is 3-D by nature, a 2.5-D model has serious limitations in space 
weather predicting, and a full 3-D numerical model is indispensable to describe realistic 
observation events. On the one hand, data-driven 3-D models can be tested and improved 
by using observation data; on the other hand, observations can be better interpreted by 
using global 3-D models. For instance, demonstrating a good match between synthetic 
and real STEREO /SECCHI images, Lugaz et al. [2009] quantitatively analyzed and well 
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explained the January 24-25 CME event by a data-driven 3-D numerical MHD model. 
Hence, two MCs' interacting in the 2.5-D model in this paper is meaningfully generalized 
to a 3-D geometry. Such model generalization and then detailed comparison with realistic 
events are out of contents in this paper and will be addressed in our near future. 

In closing, the interaction among multiple CMEs/ICMEs can be a cause of angular de- 
flection during the CME/ICME propagation. Such angular displacement, being nonlinear 
and irreversible, results in the significant responses of interplanetary dynamics and ensu- 
ing geoeffectiveness. Therefore, when successive CMEs from the solar corona are likely 
to collide with each other obliquely in the interplanetary space, the factor of potential 
deflection due to the OC should be considered for the geoeffectiveness prediction at 1 AU, 
as well as the correlation between the near-Sun and the near-Earth observations. 
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Table 1. Assortment of simulation cases of an individual MC 



Group 


Case 


(10^ km/s) 


Comment 


IPM 


Pi 


4 


Individual-Preceding MC 

{Hmc = 1) 


IFM 


Fi, F2, F3, 


4.5, 5, 6, 


Individual-Following MC 




F4, F5, Fe, 


7, 8, 9, 






F7, Fg, Fg 


10, 11, 12 
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Table 2. Assortment of simulation cases of double MCs. Note that Vmci = 400 km/s, 
dmci = 0°, tmci = hour, Hmci = 1, Hmc2 = — 1 for all thirty- four cases. 



Group 


Case 


Vmc2 

(102 km/s) 


(degree) 


(hour) 


Comment 


EID 


El, E2, E3, E4, 
E5, Ee, E7, Eg, 

Eg, Eio, Ell, E12, 
E13, Ei4, E15, E16 


6 


10 


30.1, 12.2, 44.1, 42.1, 

40.2, 37.2, 35.1, 33.1, 
31.5, 28.2, 25.1, 22.1, 
20.1, 17.1, 15.1, 10.2 


Eruption- 
Interval 
Dependence 


ADD 


Ai, A2, A3, E2,, 
A4, A5, Ag, A7, 

As, Ag, Alo 


6 


0, 3, 5, 10, 
15, 20, 25, 30, 
40, 45, 50 


12.2 


Angular- 
Difference 
Dependence 


CID 


Ci, C2, E2, C3, 
C4, C5, Cg, C7, 
Cs 


4.5, 5, 6, 7 
8, 9, 10, 11 
12 


10 


12.2 


Collision- 
Intensity 
Dependence 
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Figure 1. 
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Figure 1. Evolution of an MC2 overtaking an MCI for Case Ei, with (a)-(c) magnetic 
field magnitude B. The white solid line denotes the MC boundary. The white radial 
dashed line is along the latitude of 4.5°. Only the part of domain is adaptively plotted to 
highlight the double MCs. 
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Figure 2. 
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Figure 2. Evolution of an MC2 overtaking an MCI for Case Ei, with (d)-(f) radial 
flow speed Vr- 
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Figure 3. 
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Figure 3. Evolution of an MC2 overtaking an MCI for Case Ei, with (g)-(i) radial 
characteristic speed of fast mode c/. 
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Figure 4. Two radial profiles along Lat.= 0° and 4.5°S for Case Ei. Note that radial 
profile of B is plotted by subtracting the initial ambient value B\t=o. The solid and dashed 
lines at each profile denote the MC core and boundary, respectively. 



In-Situ Observation at LI along Lat. = 4.5°S 
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Figure 5. In-situ synthetic observations along Lat. = 4.5°S for Case Ei. Stacked from 
top to bottom are the (a) magnetic field magnitude 5, (b) elevation of magnetic field 0, 
(c) radial flow speed v^-, (d) dawn-dusk electric field VB^, and (e) Dst index. The solid 
and dashed delimiting lines denote the MC core and boundary, respectively. 



DRAFT 



May 23, 2009, 6:29pm 



DRAFT 



XIONG ET AL.: MC-MC INTERACTION AND ITS GEOEFFECTIVENESS: 2. 



X - 43 




Figure 6. 
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Figure 6. Evolution of an MC2 overtaking an MCI for Case E2, with radial flow speed 
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Figure 7. Two radial profiles along Lat.= 0° and 4.5°S for Case E2. 



DRAFT 



May 23, 2009, 6:29pm 



DRAFT 



X - 46 



XIONG ET AL.: MC-MC INTERACTION AND ITS GEOEFFECTIVENESS: 2. 



In-Situ Observation at LI along Lat. = 4.5°S 
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Figure 8. In-situ synthetic observations along Lat. = 4.5°S for Case E2. 
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Figure 9. Time dependence of MC parameters: (a) radial distance of MC core r, (b) 
latitude of MC core 9, (c) distance between botli MC cores d, (d) MCI cross section area 
Amci, (e) MC2 cross section area Amc2- In panels (a,b,d,e) the thick dashed and solid lines 
denote the preceding MCI and following MC2 in the Case E2, superimposed with the thin 
lines for the corresponding individual MC cases for contrast. In panel (c) the thick and 
thin lines represent the coupling and non-coupling conditions between two MCs. 
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Figure 10. Comparison of latitudinal distribution of Dst index among the compound- 
stream Cases El (thin sohd), E2 (thick sohd), and corresponding individual-MC cases Pi 
(dash-dotted), F3 (dashed). 
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Figure 11. Dependence of the compound-stream parameters at 1 AU on the MCl- 
MC2 eruption delay Dt {Dt = t^ci — tmci) in group EID: the (a) Sun-Earth transient 
time TT, (b) deflection angle of each MC 59, (c) total deflection angle of double MCs 
(A^ = |(56'mci| + \50mc2\), (d) distance between the two MC cores d when the MCI core 
reaches 1 AU, (e) cross section area of each MC A, (f) Dst index. The thick dashed/solid 
lines in panels (a,b,e,f) refer to the occasion of MC1/MC2 core reaching 1 AU. The thin 
dashed and solid lines in panels (a,b,e,f) denote the isolated MCI and MC2 events for 
comparison. 
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Figure 12. Dependence of the compound-stream parameters at 1 AU on the angular 
difference —D6\t=Q of the two MC eruptions in group ADD. Here —D9\t=Q = —l-D9\t=o = 

9mcl\t=0 ~ 9mc2\t=0- 
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Figure 13. Dependence of the compound-stream parameters at 1 AU on the MC2 
eruption speed Vmc2 in group CID. 
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